Beta asymmetry parameter in the decay of 114 In 
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Abstract 

The (3 asymmetry parameter A for the pure Gamow- Teller decay of In is reported. The 
low temperature nuclear orientation method was combined with a GEANT4 based simulation code 
allowing for the first time to address in detail the effects of scattering and of the magnetic field. The 
. result, A = -0.994 ± 0.010 sta t ± 0.0l0 sys t, constitutes the most accurate value for the asymmetry 

> 

parameter of a nuclear (5 transition to date. The value is in agreement with the Standard Model 

OO 

■ prediction of A = -1 and provides new limits on tensor type charged weak currents. 
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The historic measurement of the (3 asymmetry parameter A by Wu et al. |l| strikingly 
showed the violation of parity in weak interactions. Later this parameter was and still is 
extensively studied in neutron decay, e.g. to test the unitarity of the quark mixing matrix 



and to constrain the presence of right-handed currents (for reviews see [2l, l^l, 0, S] ) • In nuclear 
decays only a few measurements of the j3 asymmetry parameter were performed with the 
aim of testing symmetries of the weak interaction and/or searching for physics beyond the 
Standard Model (see e.g. j^, The most accurate result was reported in Ref. p|, i.e. 
a (purely statistical) precision of 2 % for the case of 60 Co. Still, such measurements can 
provide interesting information on possible new physics 8] . The best choice are fast (log ft 
< 5) pure Gamow- Teller (GT) transitions and the superallowed mirror transitions between 
T = 1/2 analog states 0,13], for which nuclear structure and other corrections contribute 
at most at the permille level. J n — > J n transitions between non-analog states should be 
avoided as they usually contain small isospin-forbidden Fermi contributions that originate 
from the electromagnetic interaction and modify the value of A (e.g. Ill, [12|, ll3|]). For pure 
Fermi transitions A = 0. 

Here we report a measurement of the asymmetry parameter for the 1 + — > + pure 
GT f3~ transition of 114 In (ti/2 = 72 s) with endpoint energy of 1.989 MeV and log/t = 
4.473(5) [14|. Our method combines low temperature nuclear orientation [15) with GEANT4 
simulations to address in detail, and for the first time for this method, the effects of scatter- 
ing and of the magnetic field. The result provides new information on tensor contributions 
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to the charged weak current. Other experiments are pursuing similar goals 

The electron angular distribution for (3 decay of nuclei with vector polarization J is 
written as 



W oc 



171 Pe 

1 + —b Fierz + — ■ J A 
hip hit, 



(1) 



with E e and p e the total energy and momentum of the j3 particle and m the electron rest 
mass. Expressions for the Fierz interference term bpi eTZ and the asymmetry parameter A 
are given in Refs. jj], 8]. In the Standard Model b Fierz = 0, while for a J — > J — 1 pure 
GT transition, Asm,gt((3 t ) = Tl- The observable actually determined by experiment is 
A = A/[l + {m/(E e ))b Fierz ]. 

A pure 114 In sample was obtained from the internal transition (IT) decay of 114m In (t 1 / 2 
= 49.5 d). The latter was implanted at 70 kV with a dose of 1 x 10 12 at/cm 2 into a Fe 
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foil (purity 99.99%, thickness 50 fim). The 7 spectrum showed no sign of contamination by 
another isotope. The foil was soldered at 80 °C (to prevent diffusion of In in the Fe) on the 
sample holder of a 3 He- 4 He dilution refrigerator. The latter served to polarize the nuclei by 
cooling them to millikelvin temperatures in a strong magnetic hyperfme field (in the plane 
of the foil) induced by a superconducting split-coil magnet. 

The geometry was similar to the one shown in Fig. 7 of Ref. [18|. The (3 particles 
were observed with two planar HPGe detectors [19fl with a sensitive diameter of 20 mm 
and thicknesses of 2 mm and 3 mm, respectively, placed at 0° (axial detector) and 90° 
(equatorial detector) with respect to the vertical magnetic field (i.e. orientation) axis. They 
were installed inside the 4 K radiation shield of the refrigerator at (46 ± 1) mm from the 
sample, and operating at about 10 K with an energy resolution of 4.5 keV. They were directly 
facing the sample to minimize energy loss and scattering effects. To further minimize the 
effect of scattering in the foil (a major scattering component) the plane of the foil was 
tilted 20° towards the axial detector and rotated towards the equatorial detector over 68° (a 
compromise to permit implantation of the beam). The detectors were then at angles of 20° 
and 108.5° with respect to the orientation axis. To minimize also the effect of the magnetic 
field on the j3 particle trajectories the measurements were performed in low external fields, 
i.e. B cxt = 0.046 T, 0.093 T and 0.186 T. 

The 7 rays of 114m In and of the 57 CoFe nuclear orientation thermometer [2^] were observed 
with two large volume HPGe detectors. These were placed outside the refrigerator at angles 
of 0° and 90° relative to the magnetic field axis, at (70 ± 2) mm from the sample. 



The experimental angular distribution was determined as [2l|] W{9) = N co u{9) / N warm {9) 
with 9 the angle with respect to the magnetization (orientation) direction in the Fe foil and 
Ncoid,warm{0) the count rates in a given 7 peak, or energy bin in the (3 spectrum, when the 
sample is polarized (i.e. at millikelvin temperatures; cold) or unpolarized (i.e. at about 4.2 
K; warm). Note that using a ratio of count rates reduces several systematic effects. 

The experimental angular distribution of (3 particles emitted in allowed /3 decay from 



polarized nuclei is [211 ] 



W(9) = 1 + f^-APQ cos 9 , (2) 

with v/c the initial velocity of the detected (3 particles relative to the speed of light and Q the 
solid angle correction factor. Further, P(J, fi,T, B) is the degree of nuclear polarization of 
the state with spin J and magnetic moment \x at temperature T in the polarizing magnetic 



field B = B cxt + £?hf, with Bhf = -28.7 T the internal hyperfine field which the nuclei 
feel in the Fe foil. The temperature of the sample which, via the Boltzmann distribution, 
determines the degree of nuclear polarization is obtained from the anisotropy of the 136 keV 
7 ray of the 57 CoFe thermometer. Finally, the factor / represents the fraction of nuclei that 
feel the full orienting hyperfine interaction, [iB. In the two-site model used the fraction 
(1 — /) is supposed to feel no interaction at all. The fraction / was obtained from the 
anisotropy of the 5 + — > 1 + 190 keV IT of the implanted 114m In. This was done for each field 
value as it was observed [3] that part of the saturation magnetization of the Fe foil obtained 
at 0.5 T was lost when the field was reduced to a lower value for the measurements. Since 
emission of the 190 keV 7 ray does not cause the 114 In daughter nuclei to recoil from their 
lattice position, the values obtained for the fraction / apply to 114 In as well. 

A GEANT4 based Monte-Carlo code was used to calculate the factor QcosO which 
includes the geometry of the setup, the effect of the magnetic field on the (3 particle trajec- 
tories, and (back) scattering in the source, on the sample holder and on the detector. The 
geometrical QcosO values for zero magnetic field and for the spectrum endpoint (where scat- 
tering effects are negligible) were equal to QcosO = -0.930(6) and +0.314(2) for the axial 
and equatorial detectors, respectively. 

For the energy calibration the conversion electrons from 114 In (Fig. [TJ) and the 7 rays 
from a 60 Co source were used. The Ge detectors that were used in view of the high endpoint 
energy give rise to rather significant scattering effects; the percentage of scattered events was 
found to increase from 5 % at about 1.75 MeV to 12 % at about 1.60 MeV. Therefore, the 
analysis was limited to the highest energy part of the (3 spectrum where all disturbing effects 
are minimal. The lower energy bound for the region for analysis was then set at 1.700(10) 
MeV as simulations showed that from this energy on the values of QcosO were affected by 
less than 5 % (relative) by the magnetic field, scattering, etc. (Fig. [I]). As the count rate 
for energies above 1.830(10) MeV was marginal the upper bound for analysis was set at this 
value, such that v/c = 0.9744(15). Simulations showed the values of QcosO for the region 
from 1.700 MeV to 1.830 MeV, and for the part of the spectrum above to be identical (see 
Tabled]), indicating that the conditions for scattering and magnetic field effects in the region 
used for analysis are very similar to the ones at the endpoint. The precision to which this 
holds for both detectors and the three magnetic field values, was found to correspond to 
a 0.6 % variation of A that was subsequently assigned as a systematic error related to the 
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Monte-Carlo simulations. This turned out to be the largest systematic error. 




FIG. 1: (3 anisotropy of In for the region between 1.700 and 1.830 MeV, in a field of 0.046 
T (x 2 /^ = 0.71). Upper inset: Low energy part of the j3 spectrum with 7 rays from 57 Co and 
conversion electrons from 114 In. Lower inset: High energy part of the j3 spectrum. 

Fig. [T] shows the j3 anisotropy observed in a field of 0.046 T. Results for A obtained 
from all three measurements are listed in Table [Til An overview of systematic errors is 
given in Table IHIt The energy calibration led to a negligible systematic error. The error 
related to the precision to which the geometry of the setup and the magnetic field map 
that were used in the simulations are known, was determined by repeating the analysis 
with one standard deviation modified input for the simulations, leading to systematic errors 
on A of 0.3 % and 0.4 %, respectively. The magnetic field map was calculated from the 
magnet's dimensions and properties provided by the manufacturer. The accuracy of it was 
checked by a comparison with measured field values. The difference between the calculated 
and measured values was then used as error bar. The precision to which the geometry 
of the 7 detection setup (used to determine the fraction /) was known induced a 0.13 % 
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TABLE I: Differences of the simulated values of QcosO for the energy region from 1.700 MeV to 
1.830 MeV used for analysis Canal'), and the region from 1.830 MeV up to the spectrum endpoint 



at 1.989 MeV ( 


end'), for both detectors. 




Magnetic field 




{Qcos0) anal - {Qcos0) end 




axial detector 


equatorial detector 


0.046 T 


-0.0031(42) 


-0.0032(25) 


0.093 T 


-0.0016(66) 


+0.0039(37) 



systematic uncertainty on A. The experimental uncertainty on the magnetic hyperflne 
interaction strength fiB used in the analysis translates into an uncertainty in the degree of 
nuclear orientation. This induces a systematic error of 0.16 % from fitting the (3 asymmetry, 
to obtain f'A, and an additional 0.15 % via the determination of / from the 190 keV 7 
anisotropy. Since both errors are fully correlated they are added linearly. 

Finally, analysis also accounted for the facts that 114 In, due to its short half-life, inherits 
a small part of the polarization from its precursor, 114m In, and that some nuclei may not 
yet have reached thermal equilibrium (i.e. the polarization corresponding to the sample 



temperature) in the Fe lattice when they decay 2A\, l25|, |26|]. As the fiB values for both 
isotopes that govern these effects are precisely known and,jn addition, [iB for 114 In is very 



large both effects are small and can be fully accounted for [26|, l27| . Varying then all relevant 
parameters within their error bars translated into a variation of A of at most 0.3%. 

To account for the apparent field dependence of A a systematic error of 0.4 % was added, 
determined by the average shift of the A values required to get x 2 jv = 1.0. 

Taking into account the statistical and systematic errors discussed above (Tables [IT] and 
EH) , our experimental result is A = -0.990 ± 0.010(stat) ± 0.010(syst). 

At this level of precision a small correction for the effect of recoil order terms 
necessary. Since we deal with a 1 + — > + pure GT transition only the weak magnetism, b 
in the notation of Holstein 28j], and the first class induced tensor term, d, are important. 
Further, from systematics (see e.g. H . 



3, 3(3]) we estimate bj Ac = 4.6 ± 2.7 and d/Ac = 
±(1.2 ± 1.2) with A here being the mass of the nucleus and c the GT form factor 28(. One 
then calculates a recoil correction of —0.0041(30) which leads to A = -0.994 ± 0.010(stat) 
± O.OlO(syst). 
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is 
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TABLE II: Results for A. The fraction / was obtained from the anisotropy of the 190 keV IT 7 ray 
of 114m In. Only statistical errors are listed here. The error on A is a combination of the statistics 
on the (3 anisotropy (from which the product f A was obtained) and the statistical error on the 
fraction / that was used to extract A. The error on the weighted average of A was increased by a 



factor \[y?]v = 1.55. Systematic errors are listed in Table [TTTl 



5 ext [ T ] 


fraction / 


A 


0.046 


0.734(5) 


-1.003(9) 


0.093 


0.803(8) 


-0.987(13) 


0.186 


0.874(7) 


-0.972(11) 


weighted average 




-0.990(10) 



Being in agreement with the Standard Model value 



SM 



-1 this result can now be 
used to constrain physics beyond. Our value corresponds to a lower limit of M2 = 230 
GeV/c 2 (90% C.L.) for the mass of the weak boson eigenstate W 2 that is mainly related 
to a Wr boson that couples to right-handed neutrinos J3l|. This is less stringent than 



limits from other experiments in f3 decay 
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33, 



34 



35 
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37|] . However, our result is of 



interest for tensor type charged weak currents. Assuming maximal parity violation and time 
reversal invariance for vector and axial-vector currents, one has for a J —>■ J — 1 pure GT 
f3~ transition |8j: 



A 13 ' ~ 

GT — 



aZm T fC T + C 7 
Asm H Im 

Pe 

+ J m Re [C T 



(Ee)' 



Ct 



C a 



+ Re 



\Ct\ 2 + \C t \ 2 



(3) 



2C A 

with a the fine-structure constant, Z the charge of the daughter nucleus, 7 = a/1 — (o:Z) 2 
and Cy A 7 coupling constants for vector, axial vector and tensor interactions. Primed (un- 
primed) coupling constants are for parity violating (invariant) interactions. The term with 
r ym/{E e ) comes from the Fierz interference term bpterz m Ecj- ([T]). 

With an average kinetic energy of about 1.76 MeV for the observed f3 particles, the 
sensitivity factors in Eq.Q are aZm/(p e ) = 0.084 and r ym/{E e ) = 0.209. Assuming time 
reversal invariance, i.e. Im[{CT + C' t )/Ca] = 0, as was verified at the level of about 1% in a 



TABLE III: Total error account 



Systematic effect 


Correction 
[%] 


Error 

[%] 


Energy calibration 




0.005 


Simulations of the (3 asymmetry (QcosO) 




0.60 


Geometry of the (3 detection setup 




0.3 


Magnetic field map 




0.4 


Apparent magnetic field dependence 




0.4 


Accuracy in setting the magnetic field 




0.001 


Geometry of the 7 detection setup (/) 




0.13 


Effect of error on fiB on the values of 






/ (0.15 %) and fA (0.16 %) 




0.31 


Possible incomplete relaxation of 114 In and 






polarization inherited from 114m I n 




0.3 


Induced form factors bj Ac and dj Ac 


-0.41 


0.30 


Systematics total 


-0.41 


1.0 


Statistics (see Table 111)) 




1.0 


Sum 


-0.41 


1.4 



measurement with 8 Li 38|, and to first order in Ct, Eq. ([3]) yields for time reversal invariant, 
real tensor couplings (i.e. the Fierz interference term) —0.082 < (Ct + C' t )/Ca < 0.139 (90 
%C.L.). In Fig. [2] this result is compared to limits from other experiments in nuclear (3 
decay. 

The most accurate measurement of the (3 asymmetry parameter for a nuclear (3 transition 
to date was reported. Crucial to this was the use of a GEANT based simulation code for 
this type of measurements. Our result provides limits for time reversal invariant tensor 
couplings in the weak interaction that are competitive with those from other experiments. 
A still higher sensitivity can be obtained if (3~ transitions with a lower endpoint energy are 
observed. E.g. for E e ~ 600 keV one has , -fm/E e ~ 0.4 which yields, for the same accuracy, 
two times more stringent limits. Such measurements are in progress. 
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FIG. 2: Limits (90% C.L.) on time reversal invariant tensor couplings from correlation measure- 
ments in nuclear (5 decay: the Fierz interference term bpi erz from a spectrum shape measurement 
for 114 In (with only statistical errors quoted) 
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of 6 He 
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401 ] . the f3 asymmetry parameter A in the decay 



of In (this work) and in free neutron decay [al , and the (3-v correlation coefficient a in the decays 



42[ | and Na [431 ] (see also ref. [13]). Limits for the mixed Fermi/GT transitions of the 



neutron and Na are for scalar couplings Cs = C' s = 0. 
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